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Abstract 

Ground state of the two-dimensional hard-core-boson model subjected to external magnetic 
field and quenched random chemical potential is studied numerically. In experiments, magnetic- 
field-tuned superconductor-insulator transition has already come under through investigation, 
whereas in computer simulation, only randomness-driven localization (with zero magnetic field) 
has been studied so far: The external magnetic field brings about a difficulty that the hopping 
amplitude becomes complex number (through the gauge twist), for which the quantum Monte- 
Carlo simulation fails. Here, we employ the exact diagonalization method, with which we 
demonstrate that the model does exhibit field-tuned localization transition at a certain critical 
magnetic field. At the critical point, we found that the DC conductivity is not universal, but is 
substantially larger than that of the randomness-driven localization transition at zero magnetic 
field. Our result supports recent experiment by Markovic et al. reporting an increase of the 
critical conductivity with magnetic field strengthened. 



PACS: 74.76.W Superconducting films, 71. 23. An Theories and models; localized states, 
75.40.Mg Numerical simulation studies, 68.35.Rh Phase transitions and critical phenomena, 
75.10.Nr Spin-glass and other random models. 
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1 Introduction 

Scaling argument of Abrahams, Anderson, Licciardello and Ramakrishnan [|l[ states that in two 
dimensions, infinitesimal amount of quenched randomness should drive itinerant extended states 
to localize. That is, at absolute zero temperature, the conductivity should vanish, if there exist 
any randomnesses. There are, however, some exceptions where the above description fails. 
Metal-insulator transition found in MOSFET devise [||, ^, ^, for instance, is one of recent 
hot topics arousing much attention. In the above-mentioned scaling theory, the randomness 
perturbation appears to be marginal so that some unexpected factors, for example, many-body 
interaction and external magnetic field, would possibly change the scenario. 

Suppose that there exists an attractive interaction among particles. At the ground state, the 
particles would be unstable against bose condensation so that the system would fall into a 
superconducting state. Localization transition from the superconducting state is apparently 
out of the scope of the conventional localization theory, and has been studied extensively so 
far: In experiments, the transition is observed for metallic films |[ @, ||, and Josephson- 
junction arrays ||l^, 0. In essence, ^He film adsorbed by porous media belongs to the same 



physics as well |14| , although one cannot measure its electronic conductivity. One of the 

main concerns, particularly in experiments, would be the conductivity at the localization point. 
An argument claims that localization transition occurs at a universal condition |l8| . That is, 
at the transition point, irrespective of the samples examined, the electrical conductivity stays 
w (2e)^//i with the electron charge e and the Plank constant h. A great number of experiments 
are trying to validate this issue. From theoretical view point, the criticality itself is a matter 
of interest [6[ [T^ . It is to be stressed that the localization transition occurs at the 
ground state (zero temperature). Such ground-state criticality differs significantly from that 
of the finite-temperature transitions, and belongs to a peculiar universality class: From the 
path-integral viewpoint for the partition function, d-dimensional quantum system is regarded 
as a (c? -|- l)-dimensional classical system. The system size along the imaginary-time direction 
is given by the inverse temperature, which diverges at zero temperature. Critical fluctuation 
along the imaginary-time direction (quantum fluctuation) contributes to the nature of critical 
phenomenon, giving rise to a new universality class. This extra contribution is characterized by 
the dynamical critical exponent z; a formal scaling argument on the superconductor-insulator 
transition will be found in the literatures |jl^, We stress that such arguments on criticality 
are not of mere theoretical interest, but are readily measurable in experiments; see Table ^ 

In order to tune the localization transition, either magnetic field strength or amount of 
randomness (film thickness) has to be adjusted carefully. Because the former is far more advan- 
tageous for fine tuning, numerous experiments are devoted to the magnetic-field-tuned local- 
ization transition. In Table |^, we presented a summary of experimental results. Among them, 
we would like to draw reader's attention to the recent experiment by Markovic et al. |25| p6| , 
who scanned both two tuning parameters fairly systematically. 

While in experiments, such the subtle issues have already come under through investigation, 
in computer simulation, on the contrary, it has not yet been done to incorporate the external 
magnetic field: So far, transitions driven solely by the random chemical potential (with zero 
magnetic field) have been simulated ||2^, |8[ |2^, |3l|, ^ Q . Implementation of the magnetic 
field brings about a complication such that the hopping amplitude becomes complex number. 
Because of this, quantum Monte-Carlo simulation does not work; the inefficiency would be even 
more serious than that of the so-called negative sign problem. (Note that here, the magnetic 
field is not meant to couple to the moment of spin (as in usual spin models), but it rather couples 
to the charge (hopping) through the gauge twist.) Simulation of ref. |Q is actually dealing with 
both magnetic field and randomness. In the model, however, the randomness is not of random 
chemical potential, but is of random hopping strength; the latter may lie rather out of present 
scope: Note that the chemical potential term (particle number) and the superconductivity order 
parameter (gauge) are conjugate variables. Co-existence of these mutually conjugate terms is 
the source of fascination of this physics. On the contrary, it causes the aforementioned technical 
complications. 

Here, in order to overcome the difficulty of, so to speak, 'complex-sign problem,' we employed 
the exact-diagonalization method. This method has not been used so frequently in the course of 
the studies on this issue. Perhaps, the limitation of the tractable system sizes had been worried 
over. Yet, recent extensive simulations and subsequent careful data analyses are guaranteeing 
its validity and usefulness |2^, In particular, the exact diagonalization has an advantage 
over others that the method gives dynamical response functions js^ such as the AC conductivity 
(without resorting to the maximum-entropy method, for instance). Of course, it does manage 
the complex Hamiltonian elements. 

The rest of this paper is organized as follows. In the next section, we explain our model 
Hamiltonian of the two-dimensional hard-core-boson model with magnetic field and random 
chemical potential. Details how we had chosen the gauge will be explicated there. In Section 
1^, we present numerical results. For the first time, with use of the finite-size-scaling method, 
we observe clear evidences of the onset of field-tuned localization transition. After determin- 
ing the transition point, we evaluate experimentally accessible quantities such as the critical 
conductivity and the dynamical critical exponent. The last section is devoted to summary and 
discussions. 



2 Model Hamiltonian — two-dimensional hard-core-boson 



model 

In this section, we explain our model Hamiltonian and some of its physical ingredients. The 
model we simulated is the hard-core-boson model on the square lattice subjected to external 
magnetic field and quenched random chemical potential. The Hamiltonian is given by, 



'^{^ij^ denotes the summation over all nearest neighbors on the square lattice of size L x L. 
Periodic-boundary condition is imposed as for the horizontal (z-axis) direction. The operators 
{bi} {{bl}) are the hard-core-boson annihilation (creation) operators at site i. The gauge twist 
angles {0y} are chosen such as c/jij = {iy — {L + 1)/2)ttB for the horizontal i-j pair of j right next 
to i, and = for the vertical i-j pair. Because of this choice of gauge, a particle hopping 
clockwise around a plaquette feels, as a whole, the gauge angle of irB. That is, the bosons 
are subjected to the uniform magnetic flux ttB per plaquette. Hereafter, we call B 'magnetic 
field.' (Therefore, the situation of half flux quantum per plaquette is realized at the magnetic- 
field strength B ~ 1.) Quenched random chemical potentials {Hi} distribute over the sector 
[— -s/SA, v^A] uniformly (rectangular distribution). Therefore, the strength of randomness is 
controlled by the parameter A — [iJj^Jav, where [• • -Jav denotes the random-sample average. 

Let us try to transcribe the situation realized by the Hamiltonian (^) in a manner appealing 
our intuition more vividly: Hard-core bosons are confined in a square lattice. Throughout this 
paper, we fix the particle density to be half-filled. In the absence of B and A, a statistical me- 
chanical theorem ensures that the particles are superconducting (namely, the gauge symmetry 
is broken spontaneously) at the ground state |Q; note that because of the hard-core restric- 
tion, the onset of superconductivity is a subtle issue. The random chemical potential {Hi} 
may work so as to disturb the superconductivity. The strength of the randomness is preferably 
slightly less than that of the localization-transition threshold, and so the bosons are left still 
superconducting. Because the horizontal periodic boundary condition is imposed, the lattice 
is rolled up to form a cylinder. From each end of the cylinder, two bar magnets are inserted 
so that the cylinder is subjected to magnetic field perpendicular to the cylinder surface. This 
magnetic field is supposed to destroy the superconductivity so that Anderson localization may 
occur at a certain critical magnetic field eventually. This magnetic-field-swept transition has 
not yet been simulated, whereas the randomness-driven transition has already been simulated 



extensively IgTf |2^, g9|, |30l |3l|, ^ |3§1 . For a summary of preceding simulation results at i? = 0, 
readers may consult with Introduction of ref. psf . 

It is meaningful to transform the boson Hamiltonian (|^) into the form of quantum XX 
model. With use of the mapping relations between hard-core boson and 5 = 1/2 spin ]3^ , 

bl - S+ 

b^ = Sr, (2) 
the above Hamiltonian (P is expressed in terms of the S* = 1/2 spin operators, 

^ = E i^'^-'StS- + e-'*-5-5+) +2^H.S^, (3) 

apart from a constant term. Now, we are in a position to address some notions gained from 
this spin representation: The onset of superconductivity, that is, the spontaneous breaking of 
gauge symmetry, is interpreted as the appearance of the in-plane spontaneous magnetization; 
that is, 

Mi^ = E(^»'^J + ^'^J)- (4) 



(In the above, the spin-spin correlations over the identical sites are subtracted, because they 
just yield a constant term. That subtraction improves finite-size-scaling behavior.) It is proved 
rigorously p6|| that the in-plane symmetry is spontaneous broken at the ground state of the 
two-dimensional XX model. Therefore, it is expected that at least in the vicinity oi B = 
and A = 0, our system (^) may continue to be in a superconductor phase. The perturbations 
of B and A are supposed to destroy the superconductivity. In the spin language, B introduces 
in-plane magnetic frustration at each plaquette; experts at the quantum Monte-Carlo technique 
may be convinced now why the method does not work. On the other hand, A appears to be 
random magnetic field coupling the z-component of spin. In the next section, we will study 
how MxY gets disturbed by these perturbations numerically. 



3 Numerical results 

In this section, we present numerical results. We have performed exact-diagonalization calcula- 
tion for the Hamiltonian (^) with system sizes up to L = 5. We have fixed the particle density 
to be half-filled n(~ N/L'^) ~ 0.5. For those system sizes of odd L, we carried out two sets of 
simulations for the particle numbers of = [V^ /2] and — 1 (the bracket [• • •] denotes the 

Gauss notation). Thereby, we performed interpolation calculation (least-square fitting) with 
respect to those data so as to obtain the data at n = 0.5. In the interpolation, we used the 
relation Q{n) = a{n — 0.5)^ -I- c, which means that the physics is symmetric with respect to 
n = 0.5 (owing to the particle-hole symmetry). The amount of randomness A is kept unchanged 
throughout the simulation; namely, A = 0.7. 



3.1 Field-tuned localization 

In this subsection, we will show evidences of the superconductor-insulator transition tuned by 
the magnetic field. In the analysis, the spin language (§) is used. Namely, the localization 
transition is identified as the disappearance of the XY (in-plane) magnetic order Mxy- 

In Fig. 0, we plotted the square of the in-plane magnetization per site = [(M|-y)]av/i* 
for A — 0.7 and various B. The bracket [• • -Jav denotes the random-sample average. (• • •) de- 
notes the ground-state expectation for respective random samples. The random-sample num- 
bers are 1024, 1024, 1024 and 192 for L = 2, 3, 4 and 5, respectively. As is noted above, for 
odd L, we have performed two sets of simulations. Hence, the number of random samples is 
effectively twice as many as that indicated above for odd L. From the plot, we see that the 
in-plane magnetization gets suppressed by B. However, onset of localization transition seems 
to be less transparent. Indication of the localization transition will be seized in the following 
analyses. 

In Fig. 0, we plotted the Binder parameter U |3q] for the in-plane magnetic order; that is. 
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Finite-size-scaling behavior of the Binder parameter contains information whether Mxy devel- 
ops or not: As L is enlarged, the Binder parameter should grow (be suppressed), provided that 
the order Mxy is long-ranged (short ranged). At critical point, it remains scale-invariant with 
respect to L. That is, intersection point of the Binder-parameter curves yields the location of 
the transition point. 

From the plot, we see that an intersection point locates at B w 0.1. Namely, the super- 
conductor phase survives up to the critical magnetic field B 0.1, at which the field-tuned 
localization transition takes place eventually. However, rather large corrections to finite-size 
scaling, especially for small L, prevent precise determination of the critical point. 

We carried out an alternative analysis so as to complement the above observation. We 
calculated the gauge stiffness T js^, which is defined by the elasticity with respect to the 



boundary gauge twist; 
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(6) 



where 6 denotes the gauge- twist angle through the periodic boundary, and -Eg (6) is the ground- 
state energy under the twisted boundary condition 6. In other words, T is the elasticity 
modulus against the introduction of a magnetic flux penetrating through the cylinder; note 
that our square lattice is rolled up to form a cylinder. Therefore, the stiffness should remain 
finite in the superconductor phase, while it vanishes in the Anderson-localization phase. In Fig. 
H, we plotted the scaled stiffness L^T for the same parameter range as that of Figs. |l| and|2[ 
We observe that the stiffness is suppressed by B as would be expected. In order to gain further 
information from this scaling plot, we have to consult with the formal scaling argument on T 
|l6| , 0. According to this argument, at the localization point, the stiffness should obey the 
scaling relation T ^ ]^-(d+z-2) ^ That is, the scaled stiffness ];^d+z-2Y should be scale-invariant 
at the localization transition point. Actually, from the plot, we see that the scaled stiffness with 
the particular choice of z = 2 exhibits scale invariance at critical point B 0.1; the validity of 
z = 2 will be confirmed in the analysis of Section |3.3| . Putting together the finite-size-scaling 
behaviors of U and T, we are led to the conclusion that the field-tuned localization takes place 
at i? w 0.1. In the following, we calculate experimentally accessible quantities at the point 
B = 0.1. 



3.2 Critical conductivity 

Here, we evaluate the AC conductivity at the localization point i? = 0.1 determined in the above 
subsection. As is mentioned in Introduction, an argument claims that the DC conductivity 
would be a universal value of the order « (e*)2//i with charge of one particle e* (say. Cooper 
pair). Various experiments have tried to validate this issue. However, the results are still 
remaining controversial; see Table |^. Nevertheless, it should be noted that finite conductivity 
itself is a novel feature lying out of the scope of the conventional localization theory in two 
dimensions; see Introduction. 

According to the Kubo formula, the AC conductivity is expressed in terms of the current- 
current correlation function. 



cr{iLj) — Re 



die'-*([J,(t), J,]) 



(7) 



with the current operator 
resolvent form, 
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J2i s 6"^^ '+*=^ i5a;aL ^ aj. This formula reduces to the 
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(8) 



Hence, one is forced to calculate the inverse of the Hamiltonian matrix, which is seemingly 
impossible in practice. However, Gagliano and Balseiro found that the resolvent is expanded 
into the following continued- fraction form p5| , 
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(9) 



with the coefficients generated by the following recursion relations. 



li-il 



Pf = {hm/{h^i\h~i) (/3o = o). 



(10) 



These procedures are essentially the same as those of the Lanczos tri-diagonalization. It is one 
of major advantages of the Lanczos diagonalization that one can calculate dynamical response 
functions. 

In Fig. ^, we plotted the AC conductivity at the transition point {B = 0.1 and A = 0.7); 
we had set = 0.2. The data are averaged over random samples of 1024, 1024 and 64 for 
L — 3, A and 5, respectively. Because of finite-size energy gap above the ground state, the AC 
conductivity in Fig. ^ drops in the vicinity of a; = 0. AC conductivity curve forms a peak 
beside uj = 0, and as L is enlarged, the peak position shifts toward the center uj — 0. Hence, 
for sufficiently large system sizes, the drop of cr(a; w 0) may disappear, and a Drude-like peak 
centered at a; = may emerge instead. Therefore, we regard the maximal conductivity beside 
w = as the DC conductivity for respective L. Since the DC conductivity for each L exhibits 
finite-size correction, we need to extrapolate those finite-L data to the value of thermodynamic 
limit. In Fig. ||, we plotted the finite-L DC conductivity against l/£^; the choice of this abscissa 
scale is due to the reasoning addressed in the literature Through the least-square fit, we 
obtained the DC conductivity in the thermodynamic limit such as CTc = 0.196 ± 0.001((e*)^/ft,). 
Remarkably enough, we notice that this conductivity is larger than that of the randomness- 
driven transition cfc ~ 0.135((e*)^//i) a.t B = suggesting breakdown of the universality 
of (Tc. We will discuss this point in the last section. 



3.3 Dynamical critical exponent 

In the subsection |3.l| , we have already used the relation z = 2 for the purpose of investigating 
scaling property of the gauge stiffness T (|^). In this subsection, we obtain more conclusive 
estimate of z with use of the Rieger- Young relation Q. In Fig. ||, we show the probability 
distribution of the logarithm of the first energy gap log AE over 512 random samples with L = 5 
at the transition point {B = 0.1 and A = 0.7). According to the Rieger- Young argument, the 
low-energy tail of the distribution contains an information of the dynamical critical exponent. 
In the following, we recollect their argument briefly. Because our system is disordered, the 
first energy gap AE distributes obeying a certain probability-distribution function P. We are 
considering critical phenomenon, where the physics is scale-invariant and no such particular 
energy scale that characterize the physics exists. Therefore, it is sensible to consider P(log AE') 
rather than P{AE). Low-lying excitation may be created at a peculiar part of random sample, 
where local excitation costs very little energy. Therefore, the probability may be proportional 
to the spatial volume; namely, P oc L'^AE^ with the exponent A describing the low-energy 
tail. On the other hand, the finite-size-scaling theory insists that it should be of the form P = 
P{L/^y:) — (definition of z). With use of the relation AE ~ l/^r, the distribution 

turns out to be a function of L(A£')^/^. Therefore, we obtain the relation, 

log P(log AE) = A log AE + C, (11) 

with X ^ d/z. 

In Fig. ^, we observe that the low-energy tail is, in fact, governed by that scaling relation 
with the exponent A = 1. In consequence, we obtain the relation z = d = 2, confirming the 



preliminary estimate obtained in the subsection 3.1 



4 Summary and discussions 

We have simulated the two-dimensional hard-core-boson model (|l|) subjected to the external 
magnetic field and the random chemical potential A = 0.7. Difficulty arising from the complex 
hopping amplitudes {e"^'j }, that have been preventing the application of the quantum Monte- 
Carlo method, is circumvented by the use of the exact-diagonalization scheme. For the first 
time, in computer simulation, we observed evidences of the field-tuned localization transition 
with the finite-size-scaling technique applied to the gauge order Mxy as well as the gauge 
stiffness T. After determining the location of the localization point, we carried out calculations 



of two experimentally accessible quantities, namely, the critical conductivity and the dynamical 
critical exponent. Thereby, we obtained those estimates such as — 0.196 ± 0.01((e*)^//i) 
and z = 2. First, let us make a comparison with preceding simulation results. Although no 
former results for B ^ Q are available, there have been reported a number of simulation studies 



at i? = |27|, g8|, 29, |30|, |31], |3^, M. These simulation results are settling gradually to the 
conclusions such as Uc ~ 0.135((e^//i) and z — 2. Hence, contrary to common belief, we 
are led to conclude that the critical conductivity is not universal but increases as the external 
magnetic field is strengthened. As for z, on the other hand, we found that the dynamical critical 
exponent z = 2 is retained even for B 0. 

Secondly, let us turn our attention to making a comparison with experiments. As is noted 
above, the quantities evaluated here are measurable in experiments; a summary of experimental 
results is presented in Table ^. Among various reports listed there, in particular, we would like 
to draw reader's attention to the latest exhaustive measurement by Markovic et al. p6[ , 
who scanned both magnetic- field strength and film thickness (randomness). Their experiment 
demonstrates that the critical conductivity (resistivity) does not stay universal but increases 
(decreases) gradually as B is applied. Our numerical result supports this behavior |^3|, p6[ . 
Non-universality of CTc has been arousing much attention. Meanwhile, there have been proposed 
several ideas intending to account for that |4^, 43, ^ ^ 0; for instance, it was pro- 
posed that one should take into account dissipative modes (possibly un-paired normal electrons) 
that would give rise to extra resistivity in addition to the intrinsic one. Our result shows, how- 
ever, that without resorting to such dissipative extra modes, one can account for the gradual 
increase of CTc within the framework of the pure-boson model (|l|). As for z, it appears that the 
numerical simulation result z = 2 is not in agreement with experiments. The discrepancy had 
already come out in the past studies of i? = 0. There had been proposed an attempt at altering 
the universality class |3^ ; It was claimed that an inclusion of long-range repulsion among 
particles would alter the universality class. Because in our exact-diagonalization calculation, 
tractable system sizes might not be sufhcient to implement long-range interaction, it would be 
exceedingly cumbersome to validate this issue definitely. 

In preliminary stage of our simulation, we had swept various parameter ranges so as to 
find optimal condition to observe localization transition: Note that A should not exceed a 
threshold at which particles localize already at zero magnetic field B — 0. On the other hand, 
A should not be too weak, because superconductor-insulator transition may occur at strong 
magnetic field {B ^ 0.5). For such B in close proximity to some primal fractional numbers, 
certain specific type of vortex-lattice structure becomes favored so that simulation data with 
sweeping magnetic field suffer from insystematic behaviors. Moreover, for large B, one cannot 
judge at all whether superconductivity is destroyed solely by the magnetic field or assistance 
of randomness is also significant. After sweeping various parameter ranges, we had determined 
the optimal randomness A = 0.7, and thereby, we performed extensive simulations at this 
particular condition. It would be remained for future study to gain inclusive features for the 
whole region including strong magnetic field. 
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Figure 1: Square of the in-plane magnetization (superconductivity-order parameter) (Q) is 
plotted for A = 0.7 and various B. From the plot, we see that the in-plane order gets disturbed 
gradually by the magnetic field B. 
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Figure 2: Binder parameter for the in-plane magnetic order (||) is plotted for the same pa- 
rameter range as that of Fig. 0. The intersection point of the curves gives the location of the 
localization-transition point. Hence, we found that at 5 « 0.1, localization transition takes 
place. 
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Figure 3: Scaled gauge stiffness (||) is plotted for the same parameter range as that of Fig. |l} 
The intersection point of the curves yields the critical point. Therefore, the estimate S w 0.1 
of the former analysis of Fig. || is supported. In addition, it is suggested that z — 2 holds. This 
relation will be confirmed in the subsequent analysis of Fig. 0. 




Figure 4: AC conductivity at the localization-transition point (A = 0.7 and B = 0.1) is 
plotted. Rapid drop in the vicinity of w = is due to the finite-size effect. Hence, the maximal 
conductivity beside oj = has to be regarded as the DC conductivity for each L. 
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Figure 5: DC conductivity at the localization-transition point (A = 0.7 and B = 0.1) is 
plotted against Through the least-square fit, we obtained the critical DC conductivity in 

the i — > oo limit as cTc = 0.196 ± 0.01((e*)^//i). This value turns out to be substantially larger 
than that of the randomness-driven transition ~ 0.135((e*)^//i) at B = 0. 
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Figure 6: Probability distribution of log of the first energy gap over 512 random samples with 
L = 5 at the localization-transition point (A = 0.7 and B = 0.1). The dotted line shows the 
slope of P ~ AE. From the low-energy tail of this distribution, we obtain an estimate of the 
dynamical critical exponent z = 2. 



